PROCEEDINGS 

°F SCIENCE 



O 
O 
(N 

> 
O 



> 

(N 
(N 



o 



X 



Numerical tools for the theoretical study of QCD at 
small X 



A. Sabio Vera 

Physics Department, Theory Division, CERN, CH-1211 Geneva 23, Switzerland 
E-mail: [Agustin . Sabio . Vera@cern . ch 



P. Stephens*^ 

Institute of Nuclear Physics, Polish Academy of Sciences, 
ul. Radzikowskiego 152, 31-342 Cracow, Poland 



E-mail: pstephensSannapurna . if j . edu . pi 



In this contribution we present the status of two numerical tools designed to study the small x limit 
of QCD. The first one is a Monte Carlo simulation of the BFKL evolution equation. In design of 
this approach emphasis has been placed on exploiting the linear behaviour that many variants of 
the BFKL evolution possess. This allows us to design a procedure which can be used to study 
theoretical and phenomenological aspects of different kernels. The second one is a semi-analytic 
approach to study Lipatov's effective action which describes Reggeon interactions. The study of 
the properties of this action is very complicated and we propose using a computational tool to 
handle the large amount of non-local vertices and the derivation of higher order corrections. 
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1. Introduction 

The structure of Quantum Chromodynamics is very rich. Besides the intrinsic uncertainties for 
the description of the theory in the nonperturbative region, the regime where the couphng is small 
also presents many challenges. In order to describe some processes dominated by perturbative 
dynamics it has been important to focus the efforts on different regions of phase space where loga- 
rithmic factors get enhanced and their resummation to all orders can provide valuable information. 
One example of this kind is the Regge limit of strong interacting processes. This is a regime where 
the center of mass energy, s, is very large and terms of the form as\ns are enhanced. In processes 
characterized by two large and similar transverse scales, for a range of large energies, the underly- 
ing dynamics is linear. This implies that the growth with energy of the observables can be described 
by a Unear integral equation, the so-called Balitsky-Fadin-Kuraev-Lipatov (BFKL) equation |jl|]. 
Linearity has, as a consequence, an exponential rise of the cross section with an exponent related to 
the so-called hard pomeron intercept. The LO intercept is too large when compared to experimen- 
tal data. The inclusion of the next-to-leading (NLO) order corrections ||^, where there is an extra 
power in the coupling when compared to the logs of energy, reduces the intercept to values more 
acceptable from the point of view of phenomenology. The BFKL scattering amplitudes contain a 
large amount of physical information which can be fully extracted only by numerical methods. It 
has a Poisson like structure in the number of emissions which could be useful in the study of events 
with a large multiplicity in the final state. At NLO one is sensitive to effects like the running of the 
coupling or how to best couple the Reggeized gluons to the external particles. These issues can be 
best treated in a numerical way using the methods developed in Ref. [^]. Here we will review some 
of the equations shown in these works and we will introduce a more efficient implementation, for 
alternative ideas see Ref. [Q]. The use of the tools shown here is not limited to the standard forward 
BFKL kernel but it can be extended to the nonforward case which corresponds to colour singlet 
exchange, and more exotic situations as is the case when the QCD evolution is placed in a thermal 
bath. 

Although there should be a window of large energies where BFKL effects are dominant, lin- 
earity leads to the breakdown of unitarity since the Froissart bound is not respected for asymptotic 
energies. This is connected to the well-known sentence by the late economist Kenneth E. Boulding: 
"Anyone who believes exponential growth can go on forever in a finite world is either a madman 
or an economist". Therefore at some point unitarity corrections should be taken into account. This 
is a very complicated issue from the theoretical point of view and all proposed solutions are always 
approximations. We believe an attractive approach is that proposed by L. Lipatov in Ref. A 
non-linear gauge invariant effective action was constructed which describes the interactions be- 
tween Reggeized quarks and gluons with the usual ones. The Reggeized degrees of freedom are 
represented by Wilson lines while the ^-channel emissions belong to a kinematics more general 
than the quasi-multi-Regge region needed to obtain the NLO kernel in the sense that not only a sin- 
gle pair of partons is allowed to have a fixed invariant mass but more are allowed. These "grouped" 
emissions are well separated from each other in rapidity following the usual multi-Regge kinemat- 
ics. This action incorporates unitarity in different channels and describes the interaction between 
Reggeons as well as, through quantum fluctuations around the classical solution, higher order cor- 
rections to the BFKL kernel. In this paper we propose the use of a second tool which uses the 
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semi-analytic software Effective [|T2|] to study and probe the tree-level structure of tliis effective 
action. In combination with a topology generator, the effective vertices of Reggeons coupled to 
particles can be generated automatically. 



2. Monte Carlo implementation of BFKL equations 

The first work taking full advantage of the iterative structure present in the BFKL equation to 
develop a Monte Carlo solution for the gluon Green's function was presented in Ref. More re- 
cently this method was extended to the much more complicated case of the NLO kernel in Ref. |^]. 
Here we would like to report on recent progress on a new implementation of these ideas to continue 
investigating the gluon Green's function in different physical processes. 

The LO evolution equation can be written in a very simple form as 

A/(4,4,y) ^ j d^^^e (P-a2).-^''"tV~(4+^,4,y) , (2.1) 



/ft,4,Y) = 5P)f4-4), (2.2) 



if the inital condition 



is chosen for a function related to the gluon Green's function up to a factor including the Reggeized 
gluon propagator, i.e. 



f (ka,kb,Y) = e<''^'^')^f (ka,kh,y) ■ (2.3) 



The generalization of this formula to the case of colour singlet exchange is straightforward and 
the evolution equation can be written as 



1 H 2 



dik^-X^ 



"■'In ^ ^ \ 

xe ' ^ f{ka+k,kh,q,Yy (2.4) 

The generalization for the NLO case has a very similar structure although the kernels and gluon 
Regge trajectories are more complicated. 

The numerical results are presented in fig. [l|. These correspond to the LO forward case with 
Y = 4 and ka = (50,0) GeV. Work is in progress to implement other cases at LO and NLO as well. 
In these figures the analytic solution is found from the sum of conformal spins up to 20 and the 
maximum number of emissions in the Monte Carlo implementation is 30. As can be seen there is 
still a small discrepancy as the ratio k^/kl moves away from 1. This is still under investigation. 

One of the questions we would like to answer is what is the effect of the collinear resummation 
in ^,-space proposed in Ref. on the behaviour of the Green's function when the value of the 
transverse momenta in the incoming Reggeons is very different. It is well-known that for larger 
values of the coupling this region presents an instability in terms of convergence and it would be 
important to understand how to solve this situation in a fully exclusive way by complementing the 
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Figure 1: Comparison of LO forward analytic versus MC for a fixed angles a) ^ = b) = n/A. For 
these figures we take Y = 4 and ka ~ (50,0) GeV. The source of the discrepancy, in both figures, away from 
k^/kl = 1 is still under investigation, as is the ehancement for (j) ~Q around k^/kl = 1. 



BFKL approach with all-order contributions dominant in the coUinear region. The prescription 
proposed in has a simple implementation since it suggests that the modification needed in the 
NLO kernel to introduce the coUinear improvements is to remove the term 
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in the NLO real emission kernel and replace it with 
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(2.6) 



(2.7) 



are the coUinear coefficients of the kernel in Mellin space. 

Recently a theoretical derivation of a jet definition valid at NLO has been calculated in the 
context of the NLO BFKL equation in Ref. [p^. One of the targets of our research will be to 
implement this jet definition and study jet multiplicities of fully generated events at NLO. Issues 
related to the conservation of longitudinal components in the final state emissions can be studied 
in this context. 



3. Reggeon Effective Action 

The Reggeon effective action |[TT|] contains a large quantity of physics. The problem is how to 
retrieve the interesting and possibly new physics from this action. This is still a new area of work 
and there is little experience about the best path to proceed down. We have decided to develop a 
computational tool to help us understand this action and study its behaviour. 
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The action is the standard QCD action plus a term induced by the reggeization of the gluon. 
This term can be expressed in terms of a kinetic term, describing the propagation of a Reggeon 
field, and Reggeon-gluon coupling terms 

In order to discuss these terms, we must introduce the kinematics implicit in this action. The 
action is defined in the quasi-multi-Reggeon kinematics (QMRK) approximation. In this regime 
the final states are produced in clusters, with each cluster having total momentum 2,. Cluster k is 
composed of any number of gluons, each containing momentum p~^\ Therefore we have 

/'a+Z'b = Gi + G2 + ... + e«; Q]=M],Qk = Y,pf. (3.2) 

i 

s = IPaPb = 4E^ » Si = 2QiQi+i » \ti\ = \qj\, (3.3) 

where Pa and Pb are incoming momenta. We now introduce the light-cone vector 

n'i = P^iJE ■ {n^ f =0;n+-n-=2;n^-kT=0. (3.4) 

We can now use this projection to project the Lorentz derivatives and define the inverse projected 
Lorentz derivative acting on a field 

d^ = n>id^;^A^ = —A^. (3.5) 
o± p± 

Using these definitions the kinetic term is 

^l^, = -d+R-d^R+, (3.6) 

for Reggeon field R. In the QMRK the Reggeon momentum is transverse. This is given by the 
condition d±R'^ = which is inherently included in eqn. (3^). The induced couplings between 
Reggeon fields and gluonic fields is then 

■^ind = {y [^exp (-1 |''^'a+(x0^x'+)] d^R^ix) + (- - +)| . (3.7) 

For notational simplicity we have defined 

= -iT"Al, R± = -iT"R%, (3.8) 

where T" are the SU (3) generators in the fundamental representation. The trace is over colour 
indices. This path exponentiated integral can be written in the form 

^^exp(-i|'^A+(x')^x'+) =A+£ (^A+y, (3.9) 

and can then be used to define the interaction vertices of the Reggeon field with the gluonic fields. 

We wish to be able to use the action to calculate some observables. In order to do this a 
convenient tool to have is the effective vertices which couple n Reggeon fields to m gluonic fields. 
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Figure 2: Topology of a minus Reggeon to two gluons. Zig-zag line is a Reggeon while swirly lines are 
gluons. 



This requires generating the correct topologies and then extracting from eqn. ( |3.7| ) the bare vertices 
and applying them to the topology. For this task we are using the analytic effective action analysis 



tool Effective [|12p. This tool allows us to define the Lagrangian (with the additional definitions of 
the operators) and extract the appropriate couplings from this action. We couple this to a topology 
generator to produce the tree-level effective vertices for n Reggeon and m gluons. 

We will give now an example of calculating the effective vertex which couples a Reggeon with 
minus light-cone component to two gluons. Figure ^ shows the topologies of this effective action. 
We see that this has two contributions. The first is the bare vertex which can be found in the action. 
We see, in momentum space, this is simply 



cfiRGG 



''ind 



2 ,' A+{pi)A+ip2)R-{q) , A4p,)A^ip2)R+{q) 

-iq g \ \- 



(3.10) 

Pi Pi 

where / can be the momentum of either gluon. The colour factor must also be extracted from these 
terms. Doing so we find 

-ifr'T'rAUp,)At{p2)R-{q)- 



(3.11) 



A+{pi)A+{p2)R^{q) -- 
Using the condition + = for the minus Reggeon and 

= ^{n-YA+ + ^{n+YA- +A^, 
we can sum the two contributions (choice of / in momentum) to find 

{0\^^^\A>'a{pMl{P2)R-{q)) = igq'falAn-nn 

For the other topology we have a Reggeon-gluon transition which can be found from the 
term in the induced action, and a standard QCD three-gluon vertex. The transition term simply 



(3.12) 



(3.13) 



provides a coupUng —q^^{n Y . The three-gluon couphng vertex is 

gfahc [{q-P\)vgX^ + {PY - P2)xgnv + {P2-q)^gvx\ 

Using these we find the contribution from the (T2) topology 

gfabc [{py-P2)-g^' + {-2pi -piYin-f + {2p2+PyY{ 
which can be combined with eqn. ( 3.13| ) to give the full effective vertex 



(3.14) 



(3.15) 



^ahc 



gfc 



ahc 



{pi-P2Yg^' + i-2pi-P2Yin-Y + {2p2+PiYin 



-{n-Y{n-Y 



(3.16) 
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The software package being implemented will calculate the effective vertices in the same 
manner as what is presented here. We plan to use our code to verify the other vertices found 
in [|ri|], rather than calculating them by hand. We also plan to calculate new vertices which have 
not been presented in [[TTI]. Using an automated system has the additional advantage that it may 
allow us to integrate with other codes, e.g. matrix element generators. 



4. Outlook 

We have presented the preliminary work on two computational tools. The first is a Monte 
Carlo to compute the solution to the BFKL equation for several different kernels. We also want 
this code to be able to produce fully exclusive events at NLO, but this requires still a lot of work. 
We have shown the results of our current implementation for the LO forward kernel and discussed 
the other kernels which are of interest to us. 

The second part of these proceedings introduces another computational tool designed to begin 



the analysis of the effective action proposed in Ql IQ. This action is quite complicated, but may 



include important and interesting physics. Because of the complexity to analyze the action, we 
feel a computational tool is ideally suited for this purpose. Thus we introduce the code which we 
are currently writing to analyze this action. We also present an example of the type of calculation 
which this code will perform. 

Acknowledgments: Discussions with J. Bartels and L. Lipatov are acknowledged. PS would also 
like to thank the conference organizers for their invitation and financial support. 
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